Resonant electron heating and molecular phonon cooling in single C 60 junctions 



G. Schulze 1 , K. J. Franke 1 , A. Gagliardi 2 , G. Romano 3 , C. S. Lin 2 , A. Da Rosa 2 , 

T. A. Niehaus 2 , Th. Frauenheim 2 , A. Di Carlo 3 , A. Pecchia 3 , J.L Pascual 1 
1 Inst, fur Experimentalphysik, Freie Universitdt Berlin, Arnimallee 14, H195 Berlin, Germany 
2 Bremen Center for Computational Materials Science, University of Bremen, D-28359, Germany 
3 Universita di Roma 'Tor Vergata 7 , 00133 Roma, Italy 
(Dated: March 10, 2008) 

We study heating and heat dissipation of a single Ceo molecule in the junction of a scanning 
tunneling microscope (STM) by measuring the electron current required to thermally decompose the 
fullerene cage. The power for decomposition varies with electron energy and reflects the molecular 
resonance structure. When the STM tip contacts the fullerene the molecule can sustain much larger 
currents. Transport simulations explain these effects by molecular heating due to resonant electron- 
phonon coupling and molecular cooling by vibrational decay into the tip upon contact formation. 



The paradigm of molecular electronics is the use of 
a single molecule as an electronic device pQ. This con- 
cept is sustained on the basis that a single molecule (or 
a molecular thin film) should withstand the flow of elec- 
tron current densities as large as 10 10 A/m 2 without de- 
grading. A fraction of these electrons heat the molecular 
junction through inelastic scattering with the molecule 
[2]. The temperature at the junction is a consequence of 
an equilibrium between heating due to electron flow and 
heat dissipation out of the junction. The former is domi- 
nated by the coupling of electronic molecular states with 
molecular vibrons [2j |3j 0]. The latter depends on the 
strength of the vibrational coupling between the "hot" 
molecular vibrons and the bath degrees of freedom of the 
"cold" electrodes. 

Theoretical studies predicted that current-induced 
heating in molecular junctions can be large enough to 
affect the reliability of molecular devices [2]. However, 
experimental access to this information is very limited. 
Recent studies of the thermally activated force during 
molecular detachment from a lead [5j [6] and of struc- 
tural fluctuation during attachment to it [7] reveal that 
the temperature of a molecular junction can reach sev- 
eral hundred degrees under normal working conditions, 
thus revealing that present devices work on the limit of 
practical operability [8]. Heat dissipation away from the 
junction becomes an important issue. 

In this work, we characterize the mechanisms of heat- 
ing and heat dissipation induced by the flow of current 
across a single molecule. Our approach is based on de- 
tecting the limiting electron current inducing molecular 
decomposition at varying applied source-drain bias (i.e. 
the maximum power one molecule can sustain). We use 
a low temperature scanning tunneling microscope (STM) 
to control the flow of electrons through a single Ceo 
molecule at an increasing rate until the molecule decom- 
poses. By comparing the power applied for decomposi- 
tion (Pdec) m tunneling regime and in contact with the 
STM tip we find that it depends significantly on two fac- 
tors: i) Pdec decreases when molecular resonances partic- 
ipate in the transport, evidencing that they enhance the 
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FIG. 1: (a) STM image of a 0.2 monolayer film of Ceo on 
Cu(110) (I t =l.l nA, V s =1.75 V). The intramolecular struc- 
ture (inset; It = 1.0 nA, V s =2.25 V, [31]) reveals an intrinsic 
asymmetry consistent with the adsorption orientation from 
DFT simulations (right inset), (b) Differential conductivity 
spectrum of a Ceo molecule measured using a lock-in ampli- 
fier (Rj unc t=2.2 G£l, V ac =20 mV rms). The structure in the 
spectrum is due to Ceo electronic resonances, (c) Conduc- 
tance and current vs. Z apP r plot on top of a Ceo molecule 
(V s =0.5 V). The shaded area indicates contact regime. 



heating; ii) Pd ec increases as the molecule is contacted to 
the source and drain electrodes, revealing the heat dissi- 
pation by phonon coupling to the leads. A good contact 
between the single-molecule (SM) device and the leads 
is hence an important requirement for its operation with 
large current densities. 

Our experiments are carried out in a custom-made 
ultra-high vacuum STM at a temperature of 5 K. We 
choose a Cu(110) single crystal surface because here Ceo 
adsorbs in a well-defined configuration, between 4 of the 
top most Cu atoms [9]. By annealing a sub-monolayer 
film of Ceo to 470 K we produce ordered fullerene is- 
lands with a pseudo- hexagonal structure (Fig. []Ja)), 
in which Ceo adsorbs keeping a pent agon- hexagon C-C 
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bond pointing upwards. Scanning tunneling spectroscopy 
(STS) (Fig. [TJb)) shows a clear spectroscopic finger- 
print characterized by a strong resonance at 1.5 eV above 
the Fermi level, and associated to the alignment of the 
LUMO+1 resonance (LUMO: lowest unoccupied molec- 
ular orbital). The LUMO resonance appears as a weaker 
broader peak around 0.2 eV and is partially occupied 
[TU] . A spectrum like the one shown in Fig. [TJb) is thus 
taken here as a litmus test for the integrity of the Ceo 
molecule. 

In our experiment, we approach the STM tip [11] a 
distance Z appr towards a single Ceo at constant sample 
bias (V s ) and record the current thus flowing through the 
molecule (I(Z appr )). The tunnel regime (identified here as 
the exponential regime in the I(Z appr ) plots) extends until 
the junction conductivity reaches ~ 0.03 Go (Go =77.5 
fiS). Beyond this point, a tip-Ceo contact starts to be 
formed, and the I(Z appr ) plots deviate smoothly from the 
exponential behavior [T] [12]. Ceo is very stable under 
the proximity of the STM tip. If a small positive sample 
bias (V s <0.6) is used, Ceo withstands tip contact and 
indentations of several Angstroms, holding currents close 
to 100/iA (Fig. [ijc)). 

When the value of V s is increased above 0.6 V, Ceo 
degrades during the tip approach, as we could iden- 
tify from three facts (Fig. [2|: i) a sharp discontinuity 
in the I(Z appr ) plots indicates an irreversible change in 
the molecule; ii) the height of the degraded molecule is 
typically more than 1 A lower than its neighbors; iii) 
the dI/dV s spectrum reveals that the characteristic res- 
onance structure vanishes. The precise way in which Ceo 
is degraded can not be determined exactly in our mea- 
surements. The disappearance of molecular resonances 
in the spectra hints that the most probable result is a 
rupture of the icosahedral carbon cavity. The decom- 
position is observed solely on the molecule selected for 
the indentation. Thus, we can discard that an electron- 
induced polymerization with neighboring molecules takes 
place [13] . 

A current drop as in Figs. |2|a) and[2|b) provides the 
position Zdec and current Id ec at which molecular decom- 
position occurs. Both Zd ec and Id ec depend strongly on 
the bias value V s used. In Fig. |3ja) we plot the sta- 
tistical average of Idec as a function of V s between 0.6 
V to 3.0 V. For V s = 3.0 V an electron current Idec= 7 
/iA suffices to produce the Ceo degradation. This occurs 
when the tip is still more than 1 A away from the con- 
tact position, i.e. in tunnel regime. As V s is reduced 
Idec increases gradually (with a small plateau around 1.5 
eV and 2.5 eV) and, accordingly decomposition occurs 
with the tip closer to the fullerene, but still in the tunnel 
regime (i.e. with I(Z appr ) curves as in Fig. |2ja)). For 
V s below ~1.2 V a more pronounced increase of Id ec is 
observed when approaching the tip. In this range the on- 
set of tip-Ceo contact is already detected in the I(Z appr ) 
plots (as in Fig. j2^b)) . Electron currents as high as 70 
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FIG. 2: (Color online) I(Z apP r) plots of Ceo showing molec- 
ular decomposition in (a) tunnel regime (V s =2.0 V) and (b) 
after forming a tip-Ceo contact (V s =1.0 V, acquisition time 
t=1.5s). (c) STS spectrum of a molecule before and (d) after 
an I(Zappr) plot like in panel (a). The resonance structure 
disappears, thus indicating destruction of the fullerene cage, 
(e) Island of Ceo molecules before and (f) after performing 
several I(Z apP r) events on the molecules marked in (e). Af- 
ter this, all the marked molecules show a lower height and a 
dl/dV spectrum similar to panel (d) (I t =1.0 nA, V s =2.25 V). 



H A can flow through the molecular junction for the lowest 
bias achieving degradation (V s ~ 0.6 V). In this regime 
the STM tip indents more than 2 A beyond the contact 
position. However, we can rule out a mechanical process 
of rupture because below the threshold bias of 0.6 V Ceo 
remains unaffected upon tip indentations of more than 4 
A (see Fig. [ljc)). Instead, we note that Ceo undergoes 
a thermal decomposition on metal surfaces at tempera- 
tures around 1000 K [14j [15j [16j [17] . Thus, we consider 
that the decomposition of Ceo is a current- driven ther- 
mal process, where the critical current depends crucially 
on the applied voltage. 

The two different regimes of molecular decomposition 
in tunnel and in contact as suggested by the behavior 
of Idec with V s in Fig. [3]^a) become more evident when 
we plot the power Pd ec (Pdec= Idec x V s ) applied to the 
Ceo junction for its degradation (Fig. [3^b)). In tunnel 
regime, i.e. above V s ~ 1.2 V, Pd ec amounts to ~ 20 
/iW and shows small oscillations at the LUMO+1 and 
LUMO+2 position. For lower bias Pdec increases sharply 
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FIG. 3: (Color online) (a) Statistical average of current 
reached at the point of decomposition Id ec vs. sample bias 
V s . Two sigma error bars are indicated. All the events con- 
sidered in the statistical average (180) correspond to degra- 
dation as in Fig. 2, and were done on fullerenes surrounded 
by unperturbed Ceo molecules, (b) Bias dependence of Pdec 
(squares, left scale) and distance to contact Zd ec (diamonds, 
right scale). For events in tunnel regime, Zd ec is obtained by a 
linear extrapolation. Shaded areas indicate decomposition in 
contact regime. Gray dashed lines represent normalized con- 
ductance (dI/dV s ) plots, indicating the resonance alignment. 



up to more than 50 /iW, reflecting that here the junc- 
tion can sustain larger current densities. It should be 
noted that in this regime decomposition is attained when 
the tip is in contact with the molecule. Hence, the very 
sharp change in Pd ec as compared to the oscillations in 
Pdec around higher molecular resonances suggests that 
the contact plays an important role for the decomposi- 
tion process. 

The origin of molecular dissociation in the two different 
transport regimes can be rationalized from a conceptual 
definition of molecular temperature. The electron tun- 
neling rates used to decompose the Ceo molecule (a few 
fs _1 ) are larger than typical phonon decay rates for ad- 
sorbate systems. Hence, we consider a current-induced 
vibrational heating, in which Ceo vibrations are excited 
in a non-equilibrium distribution [I8j [19] . For a certain 
set of current and bias values, the total vibrational en- 
ergy U m stored in the molecule (and hence its tempera- 
ture T m ) depends on the balance between the heat gener- 
ated by the inelastic scattering of electrons with molec- 
ular modes, and heat dissipation into the "cold" leads 
(T=5K) [2j[4j[20]. Following this picture, a larger degra- 
dation power Pdec is due to either a less effective heat 
generation or to a more effective dissipation of heat into 
the leads. 

In order to analyze and corroborate the influence of the 
resonances and the tip-molecule contact theoretically we 
have performed model transport calculations based on 
the Non-Equilibrium Green's Function (NEGF) formal- 
ism [2TJ [22] . The Ceo molecule is relaxed on a slab of 8 
layers of copper in a 4 x 5 unit cell using density func- 
tional theory calculations [23 . The resulting structure 
is in good agreement with experiments [9 . Unoccupied 
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FIG. 4: (color online) Results of theoretical simulations: (a) 
Temperature T m of the molecule vs. V s for the indicated 
distances to contact. The oscillations of T m are associated 
to the molecular resonances (dashed line), (b) The values 
Idee (diamonds, green continuous line) and Pdec (squares, blue 
broken line) are marked for a threshold temperature of 1650 
K of the Ceo molecule. The shaded area marks the contact 
regime, (c) Ratio of phonon absorption rate (A tip ) to total ab- 
sorption rate in the junction (A tot ) vs. tip-molecule distance 
(tip-molecule contact is defined as A) at 0.4 V. (d) Temper- 
ature T m at 0.4 V vs. tip-molecule distance with (blue) and 
without (red) electron-hole pair formation in the tip. 



molecular states (LUMO, LUMO+1 and LUMO+2) are 
located at 0.7 V, 1.5 V and 3.0 V above the Fermi level, 
respectively. The tip is represented by a Cu atom ad- 
sorbed on a jellium surface. 

Inelastic electron scattering in Ceo is calculated using 
the Self Consistent Born Approximation (SCBA) [25] . 
The non-equilibrium phonon population, N q (where the 
q index runs over all 174 normal modes of vibration of 
the Ceo)? is deduced from a rate equation including the 
phonon absorption (A q ) and emission probabilities (E q ) 
in the device calculated as in ref. [20] . A "model" molec- 
ular temperature, T m , can be associated to the internal 
energy, U m = J2 q ^qN q , by assuming a Bose-Einstein 
phonon population that produces the same internal en- 
ergy as the non-equilibrium population N q [24 . 

Fig. 4 (a) shows the effect of increasing the applied bias 
on the molecular temperature T m . For the experimen- 
tal values of tip-molecule distance, the applied bias leads 
to a heating of the molecule up to temperatures above 
1000 K. Such high temperatures are reached with only a 
small fraction of electrons being inelastically scattered by 
molecular vibrons (~ 10 -3 ). The heating becomes more 
effective when the LUMO resonance enters in conduc- 
tion due to resonant electron-phonon emission [4 . The 
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temperature increase is largest when all Ceo vibrational 
modes can be excited, which lie in a band of 200 meV 
width above the resonance. A similar rise in T m appears 
at bias values when higher order resonances enter into 
conduction. Contrary to this, phonon assisted tunneling 
contributes to cooling the Ceo molecule right below the 
resonances, causing in total a modulated rise in T m . 

To simulate the experimental curves of Fig. 3, we set 
a critical temperature for decomposition, Td ec , and ob- 
tain from Fig. |4ja) the corresponding set of values Vd ec 
and Zdec- From the I-V characteristics we then extract 
the corresponding Id ec - For Td ec ~ 1650 K the applied 
power is similar to the experimental values [29 . For bias 
values larger than 1.2 V the decomposition takes place 
in the tunnel regime and a step-like behavior is obtained 
for Idee- Where Id ec remains constant Pd ec increases, thus 
producing the oscillations shown in Fig. [4^b). These re- 
flect the resonant phonon cooling and heating below and 
above a resonance, respectively. Such oscillations resem- 
ble the experimental oscillations of Pd ec at the LUMO+1 
and LUMO+2, hence confirming a resonance mediated 
mechanism of molecular heating in tunneling. 

This, however, does not explain the very large increase 
in Pdec f° r sample bias below 1.2 V. Here, as in the exper- 
iment, molecular decomposition is achieved once the tip 
is in contact with Ceo- The formation of a tip-molecule 
contact enhances the dissipation of vibrations from the 
hot molecule into the cold tip. There are two possible 
mechanisms of mode quenching: vibrational decay due 
to vibron-phonon coupling and to electron-hole (e-h) pair 
excitations. Our calculations also suggest that the for- 
mer plays a small role in the cooling effect, since the 
phonon bands at the leads are much narrower than the 
200 meV wide Ceo vibrational spectrum. On metal sur- 
faces, instead, the main mechanism is the e-h excitation 
in the leads [28] . From our calculations we can estimate 
the contribution of e-h pair creation in the tip by ex- 
tracting the corresponding phonon absorption rate (A* ip ) 
[30] . Molecular mode quenching becomes only important 
in the proximity of contact formation and further inden- 
tation in the molecule as revealed by the monotonous 
increase of A* ip (Fig. 4(c)). The effect on T m can be de- 
termined by setting A* ip to zero and calculating again 
the effective temperature. Fig. |4^d) evidences that the 
temperature thus determined increases much faster when 
approaching the tip towards the molecule than in the case 
where e-h excitation in the tip is allowed. Hence, our cal- 
culations strongly support that the sharp increase in both 
Idec and Pdec result from phonon cooling upon contact 
formation. At low bias, the LUMO resonance starts to 
be removed from the conduction window and, as shown 
in Fig. Qa), the molecular temperature decreases dras- 
tically. This explains that the thermal decomposition 
cannot be achieved in the experiment below 0.6 eV. 

In summary, our results reveal that resonant tunnel- 
ing through molecular states in a single molecule device 



can generate sufficient heat to thermally decompose the 
molecular junction. For the case of a Ceo molecule on a 
Cu(110) surface, a power of only 20 fiW is sufficient. In 
order to increase the current density a molecular junction 
can sustain, it would be useful to remove molecular res- 
onances from the transport window. However, in most 
cases the molecular device properties rely crucially on 
molecular resonances, and hence cannot be engineered 
without loosing the functionality. Good contact of the 
molecule with the leads then opens the possibility for 
the single molecule device to withstand larger current 
densities. 
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